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Abstract 

This paper defines the Arrwwid number of a recursive tiling (or space-filling curve) as the 

smallest number a such that any ball Q can be covered by a tiles (or curve sections) with 

total volume 0(volume(Q)). Recursive tilings and space-filling curves with low Arrwwid 

numbers can be applied to optimise disk, memory or server access patterns when processing 

sets of points in R d . This paper presents recursive tilings and space-filling curves with 

optimal Arrwwid numbers. For d > 3, we see that regular cube tilings and space-filling 

curves cannot have optimal Arrwwid number, and we see how to construct alternatives with 

'45 better Arrwwid numbers. 

<D 

ON 1 Introduction 

^ 1.1 The problem 

Consider a set of data points in a bounded region U of M. 2 , stored on disk. A standard operation 
on such point sets is to retrieve all points that lie inside a certain query range, for example a 

'— 1 circle or a square. To prevent large delays because of disk head movements while answering 

such queries, it is desirable that the points are stored on disk in a clustered way [21 QUI 
J> [TT] [15] . Similar considerations arise when storing spatial data in certain types of distributed 

networks |19j or when scanning spatial objects to render them as a raster image; in the latter 

00 case it is desirable that the pixels that cover any particular object are scanned in a clustered 

way, so that the object does not have to be brought into cache too often |21| . For ease of 

£N) explanation, we focus on the application of clustering to storing points on disks. 

We could try to achieve a good clustering in the following way. We divide U into tiles. The 
tiles could, for example, form a regular grid of hexagons (see Figure [T]). Now we store the points 
in each tile as a contiguous block on disk. To answer a query, say with a region Q bounded by 
a circle, we compute which tiles intersect Q. For each intersecting tile, we move the disk read 
head to the position where the first point in that tile is stored, and then we retrieve all points in 
the tile by just scanning them sequentially without further delays from disk head movements. 
Since some of the tiles that intersect Q may lie partially outside Q, some of the points thus 
retrieved may be false answers: they are no answers to our query and need to be filtered out 
in post-processing. 

The approach sketched above may work pretty well, provided the following conditions are 
met: 

1. We can compute efficiently which tiles intersect Q. 

2. We can figure out efficiently where the contents of these tiles are stored on disk (for 
example, using a small index to be kept in main memory). 
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Figure 1: Sorting points into (non-)recursive tiles. Queries with small query ranges (shaded 
disk on the left) may necessitate scanning the full contents of a crowded tile. Queries with large 
ranges (shaded disk on the right) may necessitate looking up many tiles. 

3. A couple of tiles suffice to cover Q (so that we do not have to move the disk read head to 
the starting point of another tile too often); 

4. The tiles that cover Q are not much larger than Q (so that they do not contain too many 
false answers). 

These conditions are subject to a trade-off. On one extreme, we could use only one tile that 
covers all of U (making the first three conditions trivial), bring the disk head into position only 
once, and scan all data points — including lots of false answers. On the other extreme, we could 
make the tiles so small that each of them contains very few points and we get very few false 
answers, but then we may get large delays because of disk head movements when going from 
tile to tile; the first and second condition may also be harder to satisfy in this case. The best 
balance is to be found somewhere between these extremes, but where? 

If the data points are not uniformly distributed in U and we choose any fixed tile size, then 
some of the tiles may contain many more data points than others. When the query range Q 
is a small region inside such a crowded tile, we would need to retrieve all points in the tile, 
including many false answers (see Figure [T]). At the same time, there may be regions with many 
tiles with very few data points each. Queries that cover many such tiles (see Figure [T]) may 
incur an overhead for retrieving many separate tiles that is disproportionally large relative to 
the number of data points found. 

This brings us to the topic of this paper: recursive tilings. A recursive tiling is a subdivision 
of U into tiles, that are each subdivided into tiles recursively. For example, Figure [2] shows 
a recursive tiling based on trapezoids, expanded down to the level where each tile contains at 
most one data point. We store the data points in such a way that for each tile, on any level 
of recursion, the data points within that tile are stored as a contiguous sequence on the disk. 
Hopefully, if we get it right, we can now cover the region within any circle Q with a small 
set of tiles from the level of recursion where the tiles have size proportional to Q. Thus we 
would satisfy condition three and four at the same time, and because we need to retrieve only 
few tiles, we could also hope to satisfy condition one and two. Considering that moving the 
disk head once may easily cost as much time as scanning ten thousands of points from disk, 
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Figure 2: A recursive tiling into trapezoids. 



it is really important to keep the number of tiles small. Hence the topic of this study: what 
recursive tilings make it possible to cover any disk-shaped query region Q with the smallest 
possible number of tiles, while still making sure that the total size of the tiles is at most some 
constant times the size of Q? 

In this paper, we study this question for two- and three-dimensional settings, and we give 
some basic results for higher dimensions. We also study what can be achieved by controlling 
the order in which tiles are stored. Note that until now, we only required that the data points 
within each tile are stored contiguously, but we did not assume anything about the order in 
which the subtiles of any tile are stored with respect to each other. A well-chosen order may 
have the result that some of the tiles used to cover the query range are stored consecutively on 
disk, thus eliminating the need to move the disk head when going from one tile to the next. 



1.2 The Arrwwid number 

To be able to make our problem statement more precise, we define the Arrwwid number of a 
recursive tiling of a region U as follows: 

The Arrwwid number is the smallest number a such that there is a constant c such 
that any disk Q that lies entirely in U can be covered by a tiles with total area at 
most c • area(Q). 

Less formally: the Arrwwid number is the smallest number a such that any disk is covered by at 
most a relatively small tiles. We call c the cover ratio. The definition is essentially from Asano, 
Ranjan, Roos, Welzl and Widmayer [2], and we name it Arrwwid number in their honour. 

When a recursive tiling is enhanced with a recursive definition of how the subtiles within each 
tile are ordered relative to each other, the result constitutes the definition of a recursive scanning 
order or a recursive space-filling curve (this paper uses these two terms interchangeably; they 



are explained in more detail in Section 3.1). The Arrwwid number of a space- filling curve that 
fills a region U is defined as follows: 

The Arrwwid number is the smallest number a such that there is a constant c such 
that any disk Q that lies entirely in U can be covered by a curve fragments (that is, 
sets of consecutive tiles) with total area at most c • area(Q). 
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In the above definitions we could exchange disks for squares: this would only affect the 
cover ratios c but not the Arrwwid numbers a. To see this, observe that if the Arrwwid number 
is a with respect to disks, with cover ratio c, then any square Q can be covered by a disk of 
area J • area(Q) which can subsequently be covered by at most a tiles of total area at most 
c • \ ■ area(Q)- Therefore the Arrwwid number with respect to squares is at most the Arrwwid 
number with respect to disks. Furthermore, if the Arrwwid number is a with respect to squares, 
with cover ratio c, then any disk Q can be covered by a square of area - • area(Q), which can 
subsequently be covered by at most a tiles of total area at most c • - • area(Q). Therefore the 
Arrwwid number with respect to disks is at most the Arrwwid number with respect to squares. 
It follows that the Arrwwid number is the same regardless of whether disk or square ranges are 
considered. 

The above definitions naturally extend to higher dimensions by replacing disks by balls, area 
by volume, and squares by (hyper)cubes. 

1.3 Related work 

Jagadish, Kumar and Moon et al. studied how well space-filling curves succeed in keeping the 
number of fragments needed to cover a query range low P [101 EH US]. The curve quality 
measures in their work are based on the number of fragments needed to cover a query range, 
averaged over a selection of query ranges that depends on the underlying tiling. As a result 
their results can only be used to analyse space-filling curves with the same underlying tiling; 
in particular they assume a tiling that subdivides squares into four smaller squares, expanded 
down to a fixed level of recursion. This class of curves includes well-known curves such as the 
Hilbert curve [8J and Z-order, also known as Morton order or Lebesgue order |12j . However, it 
does not include, for example, Peano's curve [17], which is based on subdividing squares into 
nine squares and seems to be the curve of choice in certain applications O [21]. The work by 
Jagadish, Kumar and Moon et al. does not enable a comparison between, for example, Hilbert's 
curve and Peano's curve. 

The Arrwwid number, as defined above, does not have this limitation: it admits a comparison 
between curves with different underlying tilings. Nevertheless Asano et al., too, only studied 
curves based on the tiling with four squares per square [2]. The Arrwwid number of such a 
tiling is four. Asano et al. studied what can be achieved by controlling the order in which 
the tiles are stored: they presented an ordering scheme, the AR 2 W 2 space-filling curve, that 
guarantees that whenever four tiles are needed, at least two of them are consecutive on disk. 
Thus these four tiles can be divided into at most three sets such that the tiles within each 
set are consecutive, and thus the Arrwwid number of the AR 2 W 2 space-filling curve is three. 
Asano et al. also proved that one cannot do better: no ordering scheme of this particular tiling 
has Arrwwid number less than three. 

In practice, the effectiveness of a curve in optimising disk access time will depend on the 
cost of "jumping" to another fragment relative to the cost of scanning a false answer. None of 
the curve quality measures in the papers cited above take this into account. Bugnion et al. [4] 
studied average disk access times as a function of the cost of jumping relative to scanning. Their 
work considers scanning orders that follow square grid tilings in such a way that each tile touches 
the previous tile in the scanning order. Note that all curves mentioned so far, except Z-order, 
are like that. In such tilings, consecutive tiles may be horizontally connected (sharing a vertical 
edge), vertically connected (sharing a horizontal edge) or they may be diagonally connected 
(touching only in a vertex). Bugnion et al. analysed the performance of scanning orders under 
the assumption that they resemble random walks with a given proportion of horizontal, vertical 
and diagonal connections (and some further simplifying assumptions). Their results suggest 
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that diagonal connections (as in the AR 2 W 2 curve) are harmful for random walks, regardless 
of the cost of jumping relative to scanning. 

1.4 Results 

In this paper, we extend the scope of our knowledge on Arrwwid numbers to different tilings, 
which do not necessarily have four squares per square, and to higher dimensions. The results 
are the following: in two dimensions, no recursive tiling and no recursive space-filling curve has 
Arrwwid number better than three if the tiles are simply connected regions in the plane. This 
paper presents recursive tilings with Arrwwid number three (regardless of the ordering), as well 
as an alternative square-based space-filling curve with Arrwwid number three but without the 
diagonal connections that seem to hurt the performance of the AR 2 W 2 curve. 

In d dimensions, no recursive tiling has Arrwwid number better than d+1. We prove that 
in three dimensions, putting the tiles in a certain order will not make it possible to get below 
this bound if the tiles are convex polyhedra. There are recursive tilings with fractal-shaped tiles 
that have Arrwwid number d+1, and tilings with rectangular blocks that have Arrwwid number 
| • 2 d . In two dimensions, space-filling curves with optimal Arrwwid number can be constructed 
by defining a good order on a regular square-based tiling, but this does not generalise to three 
(or more) dimensions. Regular (Hyper)cube-based tilings have Arrwwid number 2 , and any 
space-filling curve based on such a tiling has Arrwwid number at least 2 d — 1. For d > 3 this is 
more than what can be achieved with rectangular blocks. 

2 Two-dimensional tilings 

2.1 Defining recursive tilings 

A recursive tiling of a region U (called the unit tile) in the plane is defined by a finite set of 
recursive tiling rules. Each rule specifies (i) the shape of the region to be tiled; (ii) how this 
region is tiled with a fixed number of tiles; (iii) which rules should be applied to subdivide 
each of these tiles recursively. Figure [3] shows a number of examples. Each rule is identified 
by a letter, and depicted by drawing the shape of its region, the tiles, and within the tiles, the 
letters of the rules to be applied to them; each letter is rotated and/or mirrored to reflect the 
transformations that should be applied to the subfiles of the tile. Next to each set of rules we 
see the tiling that is produced after expanding the recursion down to tiles of a few millimeters. 

Simple tilings are recursive tilings that only use one rule (Figure |3^a,b,c)). Using a term 
coined by Mandelbrot [J3], we could also call them pertilings. In contrast, composite tilings are 
recursive tilings that use multiple rules (Figure [3^d,e)). We define uniform tilings as recursive 
tilings in which all tiles have the same shape, and each rule subdivides such a shape into an 
equal number of tiles of equal size (Figure [3](a,b,d)). By square, rectangular, triangular and 
fractal tilings we mean tilings whose tiles are square, rectangular, triangular, or fractal-shaped. 
Regular recursive tilings are recursive tilings that form a fully regular grid or tiling when the 
recursion is expanded to any fixed depth (Figure [^a)). 

The size of a uniform tiling is the number of tiles in each rule. For any given (recursive or 
non-recursive) tiling, the degree of a point p in U is the maximum number of interior-disjoint 
tiles that meet in p. The vertex degree of a (recursive or non-recursive) tiling is the maximum 
degree of p over all points p £ U. 

The results in this paper apply to simple as well as composite tilings, and uniform as well 
as non-uniform tilings. We only require that the number of recursion rules is finite, that each 
rule that defines the tiling starts with a finite region, and that the transformations that may be 
applied to the rules preserve similarity. 
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Figure 3: Examples of recursive tilings, (a) A simple, uniform, regular recursive tiling, (b) A 
simple, uniform recursive tiling, (c) A simple, non-uniform recursive tiling based on the golden 
ratio 4> = |\/5+ \ - (d) Composite and uniform, (e) Composite, non-uniform. 

2.2 The Arrwwid number of a tiling: basic bounds 

Observation 1 The uniform square tiling of Figure^a) has Arrwwid number four. 

Proof: We first prove that the Arrwwid number of the uniform square tiling is at most four. 
Let r be the radius of any disk Q. Since the tiles differ in width by a factor two from one level 
of recursion to the next, there is a level of recursion where the tiles have side lengths between 
2r and 4r, and hence, area between 4r 2 and 16r 2 . The distance between any pair of horizontal 
lines of the grid formed by these tiles is at least 2r, so Q is crossed by at most one horizontal 
grid line. By the same argument, Q is crossed by at most one vertical grid line. It follows that 
at most four tiles of this grid intersect Q, so Q can always be covered by at most four tiles with 
total area at most 4 • 16r 2 = 64r 2 . Hence the Arrwwid number of the square tiling is at most 
four (with cover ratio at most 64r 2 /-7rr 2 = 64/-7r). 

It remains to show that the Arrwwid number of the uniform square tiling cannot be three 
or less. Suppose the Arrwwid number would be at most three, so there is a constant c such that 
any disk of radius r can be covered by at most three tiles with total area at most cr 2 . Consider 
a fixed value r, a tile T with area more than cr 2 , and a disk Q of radius r placed in the middle 
of this tile T, that is, centered on the spot where the four subfiles of T meet. Thus Q intersects 
each of the four subfiles of T, and the only way to cover it with less than four tiles is by covering 
it with T itself — but the area of T is larger than cr 2 . This contradicts our assumptions and 
thus proves that the Arrwwid number of the uniform square tiling is at least four. □ 
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The above arguments illustrate that there is a relation between the vertex degree of a tiling 
and its Arrwwid number: 

Observation 2 The Arrwwid number of a recursive tiling cannot be lower than its vertex degree. 

This is because any sufficiently small circle centered on a point of maximum degree would yield 
a contradiction to the assumption that the Arrwwid number is lower. 

The subdivision of the unit tile U into its subfiles creates at least one curve c in the interior 
of U. Since tile sizes decrease by at least a constant factor with each level of recursion, further 
down in the recursion c must eventually be subdivided by crossings or T-junctions with curves 
that divide the tiles on each side of c. Thus we get points on c of degree at least three. Hence 
we get the following: 

Theorem 1 Each recursive tiling of a two-dimensional region U has Arrwwid number at least 
three. 

We now ask: are there actually any recursive tilings with Arrwwid number three? In the 
next subsections, we answer this question with "yes" . 

2.3 Recursifying non-recursive tilings 

As observed above, the Arrwwid number of a tiling is at least the vertex degree of the tiling. 
Therefore, to find a tiling with Arrwwid number three, we need to find a tiling with vertex degree 
at most three. A regular hexagonal tiling has this property, but that tiling is not recursive: a 
hexagon cannot be subdivided into similar hexagons. We will now see how a recursive tiling 
can be obtained by recursively approximating a hexagon by smaller hexagons; this procedure 
will turn the boundaries of the hexagons into fractals. The technique itself is not new; however, 
I am currently not aware of any published general description such as the one given below. 

The procedure is as follows. We start with a coarse tiling which is a regular tiling of large 
hexagonal tiles. We overlay this with a fine tiling of small hexagonal tiles, such that the fine 
tiling looks the same around each large tile (Figure [4j[a)) . Now we assign each small tile to a 
large tile, in such a way that the union of the small tiles assigned to each large tile approximates 
the shape of the large tile well. To accomplish this we assign each small tile that is completely 
contained in a single large tile T to T. The remaining small tiles are assigned according to a 
tie-breaking rule that ensures that for each large tile the union of the small tiles assigned to it 
has the same size and shape. In Figure Qb) we do this as follows: when going clockwise around 
a large tile T, starting at three o'clock, we alternate between giving a small tile to a neighbour 
of T and assigning a small tile to T. All large tiles of the coarse grid are now replaced by the 
union of the small tiles assigned to them (Figure Qc)), thus adding detail to the boundaries of 
the large tiles. We now replace the small tiles by scaled copies of the large tiles (which are no 
longer hexagons), and again replace the large tiles by the unions of the small tiles assigned to 
them. This adds even more detail to the boundaries of the large tiles (Figure |4^d)). When we 
repeat this process ad infinitum, the boundaries of the tiles will converge to fractal shapes such 
that each large tile is tiled by nine scaled-down copies of itself (Figure Qe)). 

The transformation described above changed the shape of the tiles in the coarse tiling, but 
it did not change its topological structure. The vertex degree of the resulting tiling is still three, 
but unlike the original hexagonal tiling, our fractal tiling is recursive. In fact, we have produced 
a fractal tiling with Arrwwid number three (we will see how to prove such things below). 

The technique described above requires some care: it matters how well the small hexagonal 
tiles approximate the large hexagonal tiles. Figure [5] shows a recursification with seven small 
tiles per large tile, resulting in a tiling with shapes called Gosper islands [B]. Figure [6] shows 
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(a) (b) (c) 




(d) (e) 

Figure 4: How to turn a hexagonal tiling into a recursive tiling of size nine with Arrwwid 
number three. 




Figure 5: A recursive tiling of Gosper islands. 

an attempt with four small tiles per large tile. However, it fails to produce a fractal tiling 
with Arrwwid number three: the approximation of large tiles by small tiles is so bad that the 
hexagonal tiling morphs into a regular grid of rhombuses with vertex degree four. Figure [7] 
shows another failed attempt: here the result is a tiling whose tiles have disconnected interiors, 
and it contains vertices where three different tiles meet the three largest connected components 
of a fourth tile. 

Theorem 2 A recursive tiling with Gosper islands has Arrwwid number three. 

Proof: Assume that in the course tiling neighbouring vertices are at a distance u from each other. 
In the first stage of refinement, each segment e of a large tile boundary is replaced by three 
segments of length it/v7) the first and last of which make an angle of arctan g\/3 with e. The 
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Figure 6: A failed attempt to turn a hexagonal tiling into a recursive tiling of size four with 
Arrwwid number three: the Arrwwid number turns out to be four. 




Figure 7: A failed attempt to turn a hexagonal tiling into a recursive tiling of size four with 
Arrwwid number three: the Arrwwid number turns out to be four. 

new boundary of the large tile thus stays within a distance d\ = sin(arctan g\/3) -u/y/7 = f^VS 
of the old boundary. Refining the tiles i times recursively thus keeps the boundary within a 
distance di that satisfies di < f|\/3 + rfi-i/v7; for i — > oo this converges to di < 0.199 • u. Let 
s and s 1 be the radius of the smallest disk that intersects more than three large tiles in the 
original and in the recursified large tiling, respectively. We have s = u/2, and since the tile 
boundaries in the recursified tiling are within a distance of 0.199 • u from the tile boundaries in 
the original tiling, we have s' > s — 0.199 • u = 0.301 • u. 

Now consider any disk Q of radius r. There is a level of recursion in the recursified tiling 
where u, the distance between neighbouring vertices in the underlying hexagonal tiling, is 
between r/0.301 and y/7 • r/0.301. Since r < 0.301u < s', the disk Q intersects at most three 
of the tiles on this level of the recursified tiling, and each of these tiles has area §\/3 • u 2 < 
|v3 • 7r 2 /0.301 2 < 201 • r 2 (equal to the hexagon whose position they are taking). It follows 
that Q can always be covered by three tiles of total area at most 603 • r 2 < 192 • area(Q); hence 
the Arrwwid number of the Gosper tiling is three. □ 

In a similar way one could prove that the tiling of Figure |4]^e) has Arrwwid number three. 
The recursification technique can also be applied to non-hexagonal tilings. A recursified grid of 
squares with five small tiles per large tile results in a tiling of cross-like shapes (see |14j plate 
49, or [20]). Other examples with hexagons can be seen to underly the "generalised Gosper 
curves" of Akiyama et al. [I]. Figure [8] shows a recursified tiling based on a grid of squares of 
which each column is shifted by half a square's height with respect to the next column; it has 
Arrwwid number three. 
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Figure 8: A recursive tiling with Arrwwid number three based on a "shifted" grid of squares. 



2.4 An optimal rectangular tiling 

The fractal tilings are beautiful, but have one big practical disadvantage: for any given point it 
is hard to figure out in which tiles it lies. It would be quite challenging to sort any given set of 
data points into fractal tiles, so it would be hard to use such tilings to build a data structure. 
Therefore it may be interesting to look for a tiling with Arrwwid number three that uses simpler 
tiles, namely rectangles. 

In a uniform rectangular tiling, each rule defines a subdivision of a large rectangle into, say, t 
similar, smaller rectangles which are scaled by a factor \j\ft with respect to the large rectangle. 
Let a be the width/height ratio of the large rectangle; without loss of generality, assume a > 1. 

For any horizontal line through the large rectangle, the smaller rectangles along that line 
must fill the full width of the large rectangle. This implies that for fixed a and t, the equation 



cm,, 



a 



Vt 



:i) 



has at least one non-negative integer solution for the variables ti ww and rihw Similarly, consid- 
ering vertical lines we find that 

an wh + n hh = \/t (2) 

must have a non-negative integer solution. In fact, to get an Arrwwid number of three, for each 
of the variables n ww , nh w , n w h and n^h there must be a solution to the above equations where 
this variable is at least one. Otherwise the only way to subdivide the large rectangle would be 
to put the tiles in a regular rectangular grid, which would result in an Arrwwid number of four. 

Lemma 1 The size t of any uniform rectangular tiling with Arrwwid number three is a square 
number. 

Proof: For the sake of contradiction, assume that t is not a square number. Then any solution 
to Equation [2] must have n w h > 0; this holds in particular for any solution with nhh > 0. Fix 
n wh > and > such that they satisfy Equation [2] We get: 

Vi - n h h 
a = . 

Substituting a in Equation [T] yields: 



(3) 



-Vt- 



J n hh 



nhh 



Vt. 



Tlwh Tl w h n w h n w h 

Because \/i is irrational, the terms that include \ft must be equal, that is, n ww 
then gative: a contradiction. Hence t must be a square number. 



(4) 

-n hh , but 
□ 
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Figure 9: Definition of the Daun tiling, and level-two expansion. 



Lemma 2 The width/height ratio a of the rectangles in a uniform rectangular tiling of size t 
with Arrwwid number three is a rational number with numerator at most \ft and denominator 
less than \ft. 

Proof: We have a > 1, otherwise we would get a regular grid of squares with Arrwwid number 
four. Therefore the solution to Equation [l] for which n w h > 0, must have n w h < \ft and 
a = (y/i - n hh )/n wh . □ 
The above two lemmas brought an exhaustive search for increasing values of t within reach, 
trying all eligible values of a for each t. This led to the following result: 

Theorem 3 The smallest uniform rectangular tiling (with fewest tiles in the definining rules) 
with Arrwwid number three is the Daun tiling, shown in Figure^ 

Proof: We need to prove (i) the Daun tiling has Arrwwid number three, and (ii) no smaller 
uniform rectangular tiling has Arrwwid number three. The second part is tedious and boring 
but doable by hand; skeptic readers are welcome try by themselves or to request a photocopy 
of my notes. We will now prove the first part: the Daun tiling has Arrwwid number three. 

Consider a disk Q of radius r, and let k r be the level of recursion on which the tiles have 
dimensions more than 6r x 4r and at most 24r x 16r. We define a level-fc vertex, edge or tile as 
a vertex (a meeting point of three or more tiles), edge or tile that first appears in the tiling on 
the &:-th level of subdivision down from the unit rectangle. So the corners of the unit tile (the 
level-0 tile) are level-0 vertices, the remaining corners of its sixteen subfiles are level-1 vertices, 
and so forth. Consider the boundary of any tile to be labelled as follows (see Figure [l0|a) ) : 
starting from the marked corner (the one that was the lower left corner before rotation) and 
going in clockwise direction, we divide the boundary of the tile into ten sections of equal length, 
which we label B, A, B, C, D, E, D, E, F, and A, in that order. We will now prove by induction 
that there are no points of degree four in the tiling T that results from expanding the recursion 
to a depth of k r levels. Let H (k) be the following claim: 

(i) there is no point of degree four at the starting point or along the interior of an A-section 
of a level-A: tile that coincides with a B-section of another level-k tile; 

(ii) there is no point of degree four at the starting point or along the interior of a D-section 
of a level-fe tile that coincides with an E-section of another level-k tile; 

(iii) for all (x,y) G {(A, D), (A, F), (B, C), (B, E), (C, D), (C, F), (E, F)}, there is no point of 
degree four along or at any endpoint of an x-section of a level-fe tile that coincindes with 
a y-section of another level-k tile; 

(iv) subdividing a level-k tile does not create any points of degree four in its interior. 
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Figure 10: (a) Labelling the boundary of the Daun tiling, (b) Matching sections to each other. 



We first prove H(k r ). Part (iv) is trivial because level-Av tiles are not subdivided further in T. 
Part (i) is true because neither the A-section nor the B-section is subdivided further, and at 
the end of any B-section (which coincides with the starting point of the other tile's A-section) 
the tile boundary of which it is a part continues in a straight line, so that there cannot be a 
vertex of degree four there. Part (ii) and (iii) can be verified in a similar way. 

We will now prove the following for k < k r : if H(k + 1) is true, then H{k) is true. Part 
(i): As illustrated in Figure [l0|a), an A-section s at level k is subdivided into a D-section, an 
E-section, an F-section and an A-section (going clockwise around the tile) at level k + 1. A 
B-section at level k consists of a C-section, a B-section, an A-section, and a B-section (going 
counterclockwise around the tile). Thus, along an A-section of a level-fe tile that coincides with 
a B-section of another level- A; tile, a D-section is matched to a C-section at level k + 1, E to B, 
F to A, and A to B (Figure [l0|b) ) . Because H(k + 1) is true, this does not create a degree-four 
vertex at the start or anywhere along the interior of s. Part (ii) and (iii) can be verified in a 
similar way. For part (iv), it follows from H{k + 1), part (iv), that no points of degree four 
are created in the interior of level-(fe + 1) tiles, but we have to be careful about points on the 
boundary of level- (A; + 1) tiles. Figure [To| a) illustrates that the boundary sections of level- (A; + 1) 
tiles inside a level-A; tile only coincide in the combinations listed in H{k + 1), part (i), (ii) and 
(iii); therefore the interiors of these level-(A) + 1) boundary sections do not contain any points 
of degree four. All end points of these boundary sections that lie in the interior of a level-k tile 
are a starting point of a level-(/c + 1) A-section coinciding with a B-section, a starting point 
of a level-(/c + 1) D-section coinciding with an E-section, or any end point of a combination of 
level-(/c + 1) sections as listed in part (iii) of H(k + 1). Therefore, by H(k + 1), these points 
cannot have degree four either. This proves part (iv) of H{k). 

It follows by induction that H (0) holds, that is, T does not have any degree-four vertices 
in the interior of the unit rectangle. Obviously there are not any degree-four vertices on its 
boundary either. Now observe that all tiles of T consist of six squares of a regular square grid 
with line spacing more than 2r. Any disk Q of radius r intersects at most four of these squares. 
If Q does indeed intersect as much as four squares, then the number of tiles of T that intersect 
Q is equal to the degree of the grid point p shared by these four squares. As we have just 
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(a) (b) 
Figure 11: (a) Definition of the Hilbert scanning order, (b) Definition expanded by one level. 

established, the degree of any point p in T is at most three. Hence Q can be covered by at most 
three tiles of T, which have total area at most 3 • 24r • 16r = 1152r 2 = ^parea(Q). This proves 
that the Arrwwid number of the Daun tiling is three. □ 



3 Two-dimensional space-filling curves 
3.1 Defining space-filling curves 

We define a recursive scanning order as a recursive tiling in which the rules specify not only a 
subdivision into tiles and what rules to apply to the tiles, but also an order of the tiles. One 
may illustrate this by a curve that visits the tiles recursively in order. Figure [TT|a) shows an 
example: the Hilbert curve jS]. When the tiling is refined by expanding the recursion, so is 
the curve, see Figure [OJb). When defining a scanning order, the recursion within a tile can 
be rotated, mirrored and scaled as with recursive tilings. In addition it is possible to apply 
a recursive rule with reversed order; we indicate a reverse application of a recursive rule by a 



horizontal stroke above the letter that identifies the rule (see Figure 12 for an example). 

One may use a scanning order to define a mapping from the unit interval to the unit 
tile as follows. For any tile A and any x E [0,1], we define the prefix region pre(A,x) as 
the region within A that has total area x ■ area(A) and comes first in the scanning order; 
intuitively, this region can be constructed by subdividing A recursively to a sufficiently fine 
level and collecting tiles in scanning order, until tiles with a total area of x ■ area(A) have 
been collected. We can define this more precisely as follows. If x = 0, then pre(A, x) = 0. 
Otherwise let Ai,...,Ak be the subtiles of A in order, let ai,...,a,k be their areas relative to 
A, that is, aj = area(^4j)/area(^4), and let co,...,Ck be their cumulative areas relative to A, 
that is, Cj = J2)=i a j- Now let i be the largest i such that Cj_i < x. Then pve(A,x) = 
AiU...UAi^iUpre(Ai, (x— Cj_i)/oj). Let the postfix region post (.A, x) be defined as A\pre(A, x), 
and the fragment A[x, y] as post (.A, x) n pre(A, y). 

Let U be the unit tile. For fixed y, the fragment U[x, y] shrinks to a point as x approaches 
y from below; denote this point by <r^(y). Similarly, let o^(x) be the point to which U[x,y] 



shrinks when y approaches x from above; see Figure 15 for an example. Together the functions 



x 



o~+ and define a "curve" with a discontinuity (a jump) for every x such that <r+(x) ^ a 
These functions also constitute a surjective map from the unit interval to the full set of points 
within U. Thus they "fill" U: it is a space-filling curve^\ 

1 Of course real curves do not jump. When cr-\-(x) = a^(x) for all x, either function defines a proper space-filling 
curve. However, when there are jumps, our definitions only provide surjective maps from the unit interval to the 
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In this paper we identify a scanning order with the space-filling curve defined by it; we use the 
terms (scanning) order and (space- filling) curve interchangeably. The terminology introduced 
to describe tilings will also apply to curves based on those tilings. Simple curves use one rule of 
recursion; composite curves use multiple rules. Uniform or square curves are curves based on 
uniform or square tilings, respectively. 

For any fragment U[x,y], we say that (A(x) is the entry point of the fragment, while o~+(y) 
is its exit point. We say that two consecutive fragments U[x, y] and U[y,z] connect in a point p 
if p = o~^(y) = o~^(y). Since fragments of space- filling curves fill two-dimensional regions, we 
measure the size of fragments by area, not by length. 

There are other approaches to describing space-filling curves than the one chosen above. 
For example, Peano, who was the first to invent a space-filling curve, described his curve in 
an algebraic way [T7]. Other authors describe a curve by defining polygonal approximation 
of it, with a rule on how to refine each segment of the approximating polyline recursively [6]. 
Many authors specify the regions filled by fragments of the curve (like we do) together with the 
location of the entry and exit points of such fragments, but without making reversals explicit 
(for example [2 [TBI [22] ) . Since we are concerned with the use of space-filling curves as a way 
to order points in the plane, we choose a method of description that explicitly defines how to 
order the space inside a unit tile [7] . In the following, when definitions are given for curves from 
other authors, these definitions are the result of analysing the original descriptions to find a 
definition of a corresponding scanning order in our notation. 



3.2 The Arrwwid number of a space-filling curve 

Consider a space-filling curve that fills a region U. Recall that the Arrwwid number of the curve 
is the smallest number a such that there is a constant c such that any disk Q that lies entirely in 
U can be covered by a curve fragments with total area at most c-area(Q). Since every tile in the 
recursive tiling underlying a space-filling curve is by itself a fragment of the curve, the Arrwwid 
number of a curve is never more than the Arrwwid number of the underlying tiling. However, 
the Arrwwid number of a curve may be less than the Arrwwid number a of the underlying 
tiling. This can happen if, whenever a tiles are needed to cover a disk, some of these tiles are 
consecutive along the curve. To illustrate this, a space-filling curve that has a smaller Arrwwid 
number than the underlying tiling is described in the next subsection. 

3.3 Square curves with Arrwwid number three 



Theorem 4 Dekking's curve (as defined in Figure 12(a)) is a square curve with Arrwwid num- 
ber three. 



Proof: First we establish the entry and exit points of the unit square U . The locations of entry 
and exit points are determined by the transformations of the recursive rule in the different tiles. 
The curve drawn in Figure [l2|a) does not define the locations of entry and exit points; it is just 
a sketch for illustration. To prove that the entry and exit points are actually at the locations 
which the sketch may suggest, we have to study the way in which the tiles are transformed. The 
first subtile of the unit square to be visited is the one in the lower left corner. Since the recursive 

points within U, but they are not continuous, so they are not curves, and they do not match the definition of 
space-filling curves as one would find them in the mathematics literature. Z-order is an example of a scanning 
order with many jumps. Lebesgue shows how to get a proper (continuous) curve of the Z-order. In effect he 
defines pre(A, x) in another way, to include connecting pieces of curve between the tiles |12J . The same technique 
could be applied to make other curves with jumps continuous. However, for our applications it is more convenient 
to work with a representation without the connecting pieces. 
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Figure 12: (a) Definition of a scanning order following Dekking's curve [5], with marks added 
to aid in the proof of Theorem |4j (b) When p lies between two tiles (p G {a, b, c,d}), we walk to 
the closest higher-level vertex q £ {s, t} and consider p as one of the vertices of the tile that lies 
to the right, (c) Sketch of the level- two expansion of the scanning order. This curve looks quite 
different from Dekking's sketch [5], but the difference is only superficial: our sketch connects 
the centre points of the tiles, while Dekking's sketch connects the entry and exit points of the 
tiles. 

tiling and ordering rule R is applied to that tile without any rotation, reflection or reversal, the 
first subtile to be visited within it is again in the lower left corner, and this continues to hold to 
any depth of recursion. Hence, as y approaches from above, the fragment £7[0,y] shrinks to 
the lower left corner point of U, and this is therefore the entry point of U. The last subtile of U 
to be visited is the one in the lower right corner. Its exit point is what would be its entry point 
before the reversal transformation is applied; as a result of the rotation, this point is found in 
the lower right corner. Therefore the exit point of U is its lower right corner. 

Given any disk Q of radius r, consider the square grid T that is formed by subdividing U 
recursively until the squares have width at most lOr and more than 2r. Now Q intersects at 
most one horizontal line and at most one vertical line of T ■ If at most one grid line intersects Q, 
then Q lies in at most two tiles of T, with total area at most 200r 2 = ^area(Q). Otherwise 
Q is intersected by a horizontal grid line and a vertical line, which intersect in a vertex p, and 
Q intersects the four tiles of T that share vertex p. Below we prove that two of these tiles are 
consecutive in the scanning order, so that these four tiles constitute at most three space-filling 
curve fragments with total area at most 400r 2 = -^area(Q). This will prove that the space- 
filling curve has Arrwwid number at most three; together with the lower bound of Theorem [5] 
in the next subsection this proves that the Arrwwid number is in fact exactly three. 

For our proof that around any vertex p two of the tiles meeting in p are consecutive in the 
scanning order, we use the same terminology as in the proof of Theorem [3j a level-A; feature 
(vertex, edge, or tile) is one that first appears when expanding the recursion to a depth of k 
levels down from the unit tile. Now let k be the level of p, that is, p is a level-A; vertex. Now 
p is of one of two types: it either lies in the interior of a level-(A; — 1) tile, or it lies on an edge 
between two level-(/c — 1) tiles. 

If p is of the first type, one may verify in Figure |l~2"|a) (see the vertices marked by squares) 
that two of the level-A; tiles around p, say A and B, connect in p. This implies that if we expand 
the recursion further until we obtain T, the tile of T D A that touches p connects to the tile of 
T n B that touches p. 
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Figure 13: (a) Definition of the AR 2 W 2 scanning order, (b) Sketch of the level-three expansion 
(starting from rule R). 
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Figure 14: (a) Definition of the Kochel scanning order, (b) Sketch of the level-three expansion 
(starting from rule F). 



If p is of the second type, that is, it lies on an edge between two level-(A; — 1) tiles, then 
let q be the corner shared by these two tiles which is closest to p. Now consider p as a vertex 
on the boundary of the level-(A; — 1) tile that lies to the right when walking from p to q, see 
Figure [H|b). When Figure [l2|(a) depicts this tile, then p is one of the vertices marked by 
black dots (note that here we exploit the fact that the recursive rule for this curve does not use 
reflections, so that we preserve the left/right orientations for all tiles). Again, one may verify 
in Figure [l2^a) that two level-fc tiles connect in p, and hence two tiles of T connect in p. □ 

The curve analysed above was described by Dekking in the form of a recursive polygonal 
approximation [5] (he did not give a tiling-based definition or analyse its Arrwwid number). 
There are other uniform square curves with Arrwwid number three, with smaller tilings (less 
squares per square), but these are composite curves, not simple curves like Dekking's. Asano 
et al. define the AR 2 W 2 -curve [2] (see Figure |l3| and prove that it has Arrwwid number three. 
The Kochel curve (Figure 14) is another curve with Arrwwid number three, which can be proven 
to have Arrwwid number three in a way very similar to the proof for Dekking's curve. Unlike 
the j4i? 2 W^ 2 -curve and Dekking's curve, the Kochel curve has the property that consecutive 
tiles in the order always share an edge. 

In contrast to the curves mentioned above, the following square curves all have Arrwwid 
number four: Hilbert's curve [8]; Z-order, also known as Lebesgue's curve [12]; the /3^-curve [22J; 
Peano's curve |17j and all other simple uniform square curves of size nine, such as Luxburg's 
variations [13] , R-order, and Meurthe order (see Haverkort and Van Walderveen [7] for defini- 
tions in our notation). 
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3.4 A lower bound 



Above we saw that there are uniform square space-filling curves with Arrwwid number three, 
even while the underlying recursive tiling has Arrwwid number four. For the specific case of 
square space-filling curves with four tiles per recursive rule, Asano et al. proved that this is 
optimal: the Arrwwid number cannot be less than three [2]. Below we see a different proof 
technique (which we can also generalise to three dimensions later) and prove that an Arrwwid 
number of three is in fact optimal for any space-filling curve that is based on a recursive tiling 
whose tiles are topologically equivalent to disks — note that this category includes all tilings we 
have seen so far. 

Theorem 5 Any space-filling curve based on a recursive tiling with tiles that are topologically 
equivalent to disks has Arrwwid number at least three. 

Proof: Consider any subdivision of a space-filling curve, filling a unit tile U, into a set F 1 of k 
fragments U[0, xi], U[x\, X2], ?7[x/ c _ 1 , 1], such that each fragment is a simply connected region 
in the plane, that is, topologically equivalent to a disk. Let U be the unit tile's complement 
R 2 \ U, and let FbeF'U {77}. 

The boundaries of the fragments in F' form a plane graph Q with face set F. As the edge set 
E of Q we take the maximal curves that form a boundary between two faces of F. The vertex 
set V of Q is the set of points where three or more faces of F meet. By degr(v) we denote the 
number of edges of E that are incident on the vertex v. By ends(u) we denote the number of 
fragments of F' that have v as their entry point plus the number of fragments of F' that have 
v as their exit point. 

Let V be the vertices of Q that are not on the outer face U. I claim that if the space-filling 
curve has Arrwwid number less than three, we must have degr(-u) — ^ends^) < 2 for every 
vertex v G V . To prove this, let's assume, for the sake of contradiction, that there is a vertex 
v € V' with degr(u) — ^ends^) > 2. Let F\, ...,F m be the faces of F that meet at v (where 
m = degr(u)), in the order in which they appear in the scanning order; each face F{ constitutes 
a fragment U[fi, f'A. For 1 < i < m, let Bi (bridge i) be the smallest fragment C/[6j,6^] (not 
necessarily a fragment in F') that starts from v in F{ and ends at v in i^+i, more precisely: B{ 
is the smallest fragment U[bi, 6'J such that bi € f-], b\ 6 [fi+i, and cr^(bi) = ^(b'A = v 



(see Figure 15 for an example). Note that Bi degenerates to the point v if and only if Fi and 
Fi + \ connect at v, contributing two to the entry and exit points counted in ends(f). Therefore 
the number of degenerate bridges at v is at most |ends(t>). Since ^ends(v) < degr(v) — 2 and 
there are degr(t;) — 1 bridges in total, there are at least two non-degenerate bridges at v. By 
assumption the space-filling curve has Arrwwid number less than three, so there must be a 
constant c such that any circular query range Q can be covered by at most two fragments (not 
necessarily from F') with total area at most c • area(Q)- Now consider a circle Q centered at v 
with area less than 1/c times the area of the smallest non-degenerate bridge at v. To cover Q 
with only two fragments, the fragments must cover Q and all but one of the bridges B\, B m -\. 
Thus at least one non-degenerate bridge is covered, and by the definition of Q, this bridge has 
area more than c • area(Q). However, this contradicts the definition of c. Therefore, if the 
space-filling curve has Arrwwid number less than three, we must have degr(f ) — ^ends(t> ) < 2 
for every vertex v G V . 

Summing the above claim over all vertices in V we get: 

degr(w) < 2\V'\ + 1 ^ ends(v) < 2\V\ + 1 ^ ends(w). (5) 

v&v vev vev 

Let E' be the set of edges with both end points in V , and let E" be the remaining edges (those 
with at least one end point on the outer face). We have YlveV degr(v) > 2\E'\ = 2\E\ — 2\E"\ 
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Figure 15: Illustrating the proof of Theorem [5j On top we see a unit square divided into sixteen 
tiles, and below we see the corresponding unit interval that represents the space-filling curve, 
divided into sixteen fragments. The four tiles meeting in v are F\,F%, F% and F4, and the bridges 
around v are B\ (dark shaded area), B2 (light shaded area), and B3 (just the point v). Any 
disk Q centered on v can be covered with three small fragments (two squares, and one rectangle 
consisting of two squares). To do with fewer fragments, that is, only two, these fragments need 
to cover at least one of the bridges B\ and B2 ■ 

and ^2 v ^yi ends(v) < 2\F\ — 2 (since each face except the outer face has one entry point and 
one exit point). Thus Equation [5] implies: 

2\E\ - 2\E"\ < 2\V\ + \F\ - 1. 

Moving 2\E"\ and F to the other side and subtracting Euler's formula \E\ — \F\ = \V\—2 twice 
we get: 

\F\ < 2\E"\ + 3. (6) 

Now take any space-filling curve based on a recursive tiling whose tiles are simply connected 
regions (topological disks). Consider the subdivision of the unit tile U into tiles down to a 
level of recursion on which there is a tile T that does not touch the boundary of U. This level 
must exist, since the maximum diameter of the tiles decreases geometrically with the level of 
recursion, so at some point the tiles become so small that the tile that covers the centre point 
of U does not touch the boundary of U. Now subdivide T into smaller tiles recursively while 
keeping the tiles outside T as they are. Thus we keep E" fixed while we continue to increase 
\F\, eventually getting \F\ > 2\E"\ + 3. This contradicts Equation [6j so the space-filling curve 
cannot have an Arrwwid number less than three. □ 
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Figure 16: Approximating large cubes in a tiling with vertex degree four by 125 small cubes 
each, taken from the same tiling scaled by a factor 1/5. 



A direct consequence of the above is that when a space-filling curve is based on a recursive 
tiling with simply connected tiles and Arrwwid number three, the Arrwwid number of the space- 
filling curve is always exactly three as well, regardless of the order in which the curve traverses 
the tiles; the "curve" does not even need to be continuous. 



4 Three-dimensional tilings 

The definitions of tilings and Arrwwid numbers generalise naturally to higher dimensions. The 
Arrwwid number of a three-dimensional recursive tiling of a unit tile U is the smallest number 
a such that there is a constant c such that any ball Q that lies entirely in U can be covered 
by a tiles with total volume at most c • volume(Q). Theorem [T] generalises naturally to three 
dimensions: 

Theorem 6 Each recursive tiling of a three-dimensional region U has Arrwwid number at least 
four (proof omitted). 



4.1 Optimal recursified tilings 

This subsection describes how to construct a recursified three-dimensional tiling with optimal 
Arrwwid number, that is, four. The construction is based on a tiling with cubes that are shifted 



with respect to each other (Figure 16). The coarse tiling consists of horizontal layers that are 
one cube high; each layer is shifted to the right and to the front with respect to the layer below 
over 1/3 of a cube's width. Each layer consists of columns that are one cube wide; each column 
is shifted to the back with respect to the column to the left over 1/3 of a cube's width. The 
fine tiling is equal to the coarse tiling scaled by a factor 1/5 (with the centre point of a large 
cube as its fixed point); each small tile is assigned to the large tile with which it has the largest 
overlap. 

Theorem 7 There is a recursified three-dimensional tiling with Arrwwid number four. 

Proof: The proof follows the approach of Theorem [2] The recursified tiling is constructed from 
the coarse and fine tilings of shifted cubes as described above. Let w be the width of a cube 
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Figure 17: Lifting a Daun tiling into three dimensions. 



in the coarse tiling, and thus, w' = w/5 is the width of a cube in the fine tiling. In the coarse 
tiling, the smallest ball that intersects more than four tiles has radius w/6. When we replace a 
large tile by the union of 125 small tiles, the boundary of the large tile stays within a distance of 
^\/2-w' = =jij\/2"it> from its original location. In recursion with scale factor 1/5, the movement 
of the boundary adds up to at most | • j^V? ■ w = j^V^ ■ w. Thus the smallest ball that 
intersects more than four large tiles will still have radius at least (g — j2\/2)w > ^w- The 
proof can now be completed as in Theorem [2] □ 
Another way to get this result is based on a regular tiling with truncated octahedra. In such 
a tiling all vertices are incident on four tiles. The coarse tiling contains the truncated octahedron 
whose vertices have coordinates (0, ±5, ±10), and all permutations of these coordinates. It is 
the intersection of an axis-parallel cube with diagonal (—10,-10,-10) — (10,10,10), and an 
octahedron with vertices (±15,0,0), (0,±15,0), and (0,0, ±15). Further truncated octahedra 
are placed at translations (20fc, 20/, 20m) and (20fc + 10, 20/ + 10, 20m + 10), for all k,l,m£ Z. 
The fine tiling is the same as the coarse tiling scaled by a factor 1 /5. To the large tile centered at 
the origin we assign 125 small tiles, namely those with centre points at the following coordinates 
or permutations thereof: (0,0,0), (0,0, ±4), (±2, ±2, ±2), (0,±4,±4), (0,0, ±8), (±2, ±2, ±6), 
(±4, ±4, ±4), (±2, ±6, ±6), (±0,±4,±8), (2,2,10), (-2,-2,10), and (-2,2,-10). To prove 
that the resulting tiling has Arrwwid number four, we need the following. In the coarse tiling, 
the smallest ball that intersects more than four tiles has radius §\/2. When we replace a large 
tile by a union of small tiles, the boundary of the large tile stays within a distance of 2 from 
the original boundary. In recursion with scale factor 1/5, the movement of the boundary adds 
up to at most | • 2 = I. Thus the smallest ball that intersects more than four large tiles will 
still have radius at least ^(V2 — 1). The proof can now be completed as in Theorem [2J 

4.2 Rectangular tilings 

The following is a straightforward generalisation of Observation [T] 

Observation 3 Any uniform cube tiling has Arrwwid number eight. 

With rectangular tiles, that is, with axis-parallel boxes, one can do better. The three- 
dimensional rectangular tiling with lowest Arrwwid number found so far, is the "lifted Daun 



tiling" shown in Figure 17 It is obtained by adding a third dimension to the two-dimensional 
Daun tiling from Figure [9} The unit tile of the three-dimensional tiling is a rectangular block 
with width-to-depth ratio 3/2 and arbitrary height. It is divided into 64 subfiles, organised 
into 4 equal layers of 16 tiles each. Each layer, seen from above, shows a Daun tiling. The 



20 



tiles are rotated around a vertical axis only in the same way as in the Daun tiling. Each vertex 
in the resulting tiling is on the boundary between two layers and is adjacent to three tiles in 
each layer (because the Daun tiling has vertex degree three), therefore the vertex degree of 
the three-dimensional tiling is six. From here one can prove (similar to Theorem [3| that the 
three-dimensional "lifted Daun tiling" has Arrwwid number six. 

Theorem 8 There is a three-dimensional rectangular recursive tiling with Arrwwid number six. 

We do not know if it is possible to get an Arrwwid number of four or five with a rectangular 
tiling in three dimensions. As in the two-dimensional case, we could derive some properties that 
the aspect ratios and the numbers of tiles should have in order to allow an Arrwwid number 
of four. Nevertheless the search space is still huge. So far I only managed to search all three- 
dimensional uniform rectangular tilings with less than 27 tiles; no recursive tiling with Arrwwid 
number four was found. 

5 Three-dimensional space-filling curves 
5.1 Rectangular space-filling curves 

Theorem 9 Any space-filling curve based on a uniform cube tiling has Arrwwid number at least 
seven. 

Proof: Consider the regular cube tiling T obtained by subdividing a unit cube U recursively 
into smaller cubes to a certain depth of recursion. Let C be the set of tiles obtained. By tiles(t;) 
we denote the number of tiles in C that are incident on the vertex v. By ends(-u) we denote the 
number of tiles in C that have v as their entry point plus the number of tiles in C that have v 
as their exit point. 

Let V' be the set of vertices of T that do not lie on the boundary of U. If the space-filling 
curve has Arrwwid number less than seven, we must have tiles(u) — ^ends(v) < 6 for every vertex 
v G V' (the proof of this claim is a straightforward adaptation of the proof of Theorem [5]). Note 
that tiles(v) = 8 for each v G V, so in fact we must have ^ends(u) > 2 for each v G V . Since 
the total number of entry and exit points at vertices is at most 2|C|, this leads to: 

ici > i Yl ends W ^ E r nds(w) - 2|n (7) 

However, in a uniform cube tiling we have \V'\ = (ICI 1 / 3 - l) 3 , which is more than \\C\ for 
a sufficiently deep level of recursion. This contradicts Equation [7j so the space-filling curve 
cannot have an Arrwwid number less than seven. □ 

Theorem 10 Any space-filling curve based on the lifted Daun tiling of Section \4-S\ has Arrwwid 
number six. 

Proof: The proof is similar to that of Theorem [9} In this case a contradiction is derived from 
tiles(f) — ^ends(w) < 5 and tiles(v) = 6 for each v G V, which gives |ends(t>) > 1 for each 
v G V . Summing up over all v G V we get \C\ > \V'\. Now consider the tiling obtained by 
expanding one level of recursion of the lifted Daun tiling. Whenever we subdivide a tile in this 
tiling, we replace a tile by 64 smaller tiles, and add a number of new vertices: to start with, on 
each of the three boundaries between the four layers of new tiles, there are 30 new vertices (see 
Figure [9| , at least 24 of which lie in the interior of the unit tile U (since at least one short side 
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Figure 18: For a vertex v on a polyhedron t, we have that Angle(v, t) and Turn(i>, t) can be 
covered by the triangles between them on the unit sphere of directions. 

and at least one long side of the tile lies in the interior of U). Thus we increase the number of 
tiles by 63 while increasing the number of vertices in the interior of U by at least 72. Therefore, 
when we subdivide enough tiles, we get \V'\ > \C\, which contradicts the above. Therefore the 
space-filling curve cannot have Arrwwid number less than six. □ 

5.2 A lower bound for convex space-filling curves 

To prove a lower bound on the Arrwwid number of three-dimensional space-filling curves, we 
follow the same general approach as in two dimensions. The proof in the two-dimensional case 
essentially has two main ingredients: (1) In any subdivision of a two-dimensional unit tile into 
simply connected tiles, the total number of vertex-tile incidences is roughly twice the number 
of edges; therefore Euler's formula gives us a relation A between the number of vertices, the 
number of tiles, and the total number of vertex-tile incidences. (2) For the Arrwwid number 
to be lower than three, there must be another relation B between the number of vertices, the 
number of tiles, and the total number of vertex-tile incidences. It turns out A conflicts with B 
for large enough tilings, hence an Arrwwid number lower than three is impossible. 

Unfortunately ingredient (1) does not easily generalise to three dimensions, because be- 
sides the number of edges, another unknown enters the equations, namely the number of two- 
dimensional facets between the tiles. Therefore, to be able to complete our proof, we need 
another way to establish a relation between the number of vertices, the number of tiles, and the 
number of vertex-tile incidences in a three-dimensional tiling. We will now see how this can be 
done when the tiles are convex polyhedra. 

Consider a convex polyhedron t, and let v be a point on the boundary of t. In the following, 
we define the size of a set of vectors S as the area of the projection of the vectors in S on a 
unit sphere. Let angle(w,t) be the size of the set Angle(v,t), which consists of the vectors that 
point from v into t. Let turn(f , t) be the size of the set Turn(u, t), which consists of the vectors 
that point from v away from t at an angle of at least 7r/2 with the boundary of t. Figure [18^ a) 
illustrates the case in which v is a vertex of t. Note that if v is a point in the interior of an edge 
of t, we have angle (u , t) G (0, 2ir) and turn(t>, t) = 0; if v is a point in the interior of a facet of t, 
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we have angle(u, t) = 2tt and turn(w, t) = 0. We now have the following: 

Lemma 3 turn(i>, t) + angle (u , t) < 2tt. Equality holds if and only if v is a point in the interior 
of a facet of t, in which case angle (u , t) = 2tt and turn(t>, t) = 0. 

Proof: Let ei, e m be the edges incident to v. For any edge ej, let h(&i) be the plane through 
v that is orthogonal to e«, and let H(ei) be the halfspace bounded by /i(e$) that does not 
contain ej. Observe that the vectors that point from v away from t at an angle of at least tt/2 
are exactly those that point into the intersection of the halfspaces H(e\), ...,H(e m ). Thus, on 
the unit sphere of directions, we find that Angle(v, t) is represented as a spherical polygon whose 
vertices correspond to the directions of the edges e±, ...,e m with respect to v, and Turn(u,t) is 
represented as a spherical polygon whose edges are segments of h(ei), h(e m ), see Figure |l8|h). 
The relation between a and h(&j) as explained above implies that the remaining part of the 
unit sphere, between Angle(v,t) and Turn(w,t), can be triangulated by arcs of length ir/2 that 
make right angles with the adjacent edges of Angle(t>,£) and Turn(u, t). 

Now consider Turn(t>,t) and the adjacent triangles, see Figure |l8fc). Since Turn(v,t) is 
smaller than a hemisphere, the triangles would cover Turn(w, t) completely if folded onto it, see 



Figure 18 d). Similarly, the triangles adjacent to Angle (u, t) would cover Angle(v, t) completely. 
Hence the total area of the triangles is at least angle(t> , t) + turn(u, t), so angle(u, t) + turn(u, t) 
is at most 2ir, the area of half a unit sphere. 

Equality holds only if the triangles adjacent to Turn(t>, t) cover Turn(u, t) exactly, and those 
adjacent to Angle(t>,i) cover Angle(u,i) exactly. This only happens when both Angle(u,£) and 
Turn(u, t) are either empty, or a full hemisphere. Angle(v , t) is never empty, but it is possible 
that angle(u, t) = 2tt and tum(v, t) = 0, namely if v is a point in the interior of a facet of t. □ 

Consider a subdivision of a bounded convex polyhedral unit tile U into a set C of convex 
polyhedral tiles, whose shapes come from a fixed set S. Let V be the set of vertices of the 
subdivision, and let V' C V be the vertices that lie in the interior of U. We can now prove: 

Lemma 4 \V'\ + (1 + a)|C| < h X^ueV tiles(u), where tiles(f) is the number of tiles with v on 
their boundaries, and a > is a constant that depends on S. 

Proof: For ease of notation, define angle(u,t) = turn(w,t) = when v does not lie on the 
boundary of t. Observe that for every vertex v £ V' the sets Angle(f , t) of all incident tiles t 
together cover the full sphere of directions, so we have ^ tgC - angle (u, t) = 4ir. Also, for every 
tile t G C the sets Turn(u , t) of all vertices v on t together cover the full sphere of directions, so 
we have X^eV t urn ( u > t) = 4tt. With these observations we get: 

\V'\ + \C\ = 

h E«eV ^tec angle(u, t) + ± J2teC J2vev turn(u, t) < 
h T^vev Etec (angle(v, t) + turn(w, t)) . 

Since every tile t G C is a bounded convex polyhedron with a shape from a fixed set S, it 
must have at least one vertex v with tum(v,t) > 0. Hence, by Lemma pjl for this vertex v we 
have angle(u,i) + turn(f,t) < 2ir — (3, where /3 > is a constant that depends on S. For other 
vertices v on the boundary of t we have angle(f ,t) + turn(w, t) < 2ir (by Lemma [3]), and for 
vertices v that are not on the boundary of t we have angle(-u, t) + turn(v, t) = (by definition). 
Therefore the above equation gives: 



V v&v / " \ vev 



Setting a = /3/4-7T proves the lemma. □ 
We are now ready to prove the lower bound we are after: 
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Theorem 11 Any space-filling curve based on a recursive tiling with convex tiles in three di- 
mensions has Arrwwid number at least four. 

Proof: Consider a subdivision T of the unit tile U into a set C of convex polyhedral tiles, 
obtained by applying the recursive rules that define the tiling that underlies the space-filling 
curve. Let V be the set of vertices of T, and let V' C V be the vertices that are not on the 
boundary of U. Let tiles(u) be the number of tiles in C that are incident on the vertex v. Let 
ends(w) be the number of tiles in C that have v as their entry point plus the number of tiles in 
C that have v as their exit point. 

If the space-filling curve has Arrwwid number less than four, we must have tiles(u) — 
^ends(i>) < 3 for every vertex v G V (again, the proof of this claim is a straightforward 
adaptation of the proof of Theorem [5]). Summing over all vertices in V' gives: 

tiles(» < 3\V'\ + 1 ^2 endsO). 

Let V" be the set of vertices of T on the boundary of U, that is, V" = V \ V . Now Lemma [4] 
gives \ V'\ + (1 + a)|C| — \ YlveV" tiles(u) < \ YlveV tiles(u) for some fixed constant a. Because 
all tiles are convex polyhedra, all vertices v £ V have tiles(u) > 4, and therefore 2\V'\ < 
\ YlveV tiles(u). Note that the total number of entry and exit points at vertices is at most 
2|C|. Therefore the above equation gives us: 

3\V'\ + (1 + a)\C\ - g tiles (^) < 2 tiles W ^ 3 \ V '\ + lYl ends ( y ) ^ 3 \ V '\ + \°\- 

Therefore we must have: 

a \°\ < \ E tiles ( v )- ( 8 ) 

v&V" 

From here we follow the same approach as in the proof of Theorem [5| we construct T by 
subdividing the unit tile U into tiles recursively until we get a tile T that does not touch 
the boundary of U. Then we subdivide T further, keeping the vertices V" and their incident 
tiles fixed, while continuing to increase |C|, eventually getting a\C\ > \ YlveV" tiles(u). This 
contradicts Equation [8j so the space-filling curve cannot have Arrwwid number less than four. 

□ 

6 More than three dimensions 

Many of the results for three dimensions, but not all, generalise to higher dimensions in a 
straightforward way. These results are given here without writing out the proofs: 

Theorem 12 Each recursive tiling of a d- dimensional region U has Arrwwid number at least 
d + 1. 

Theorem 13 There is a recursified d-dimensional tiling with Arrwwid number d+1. 
(The tiling is constructed from a tiling of shifted hypercubes as in Section [4. 1| ) 

Theorem 14 Any uniform hypercube tiling has Arrwwid number 2 d , and any space-filling curve 
based on it has Arrwwid number 2 d — 1 or 2 d . 

Theorem 15 There is a rectangular recursive tiling in d dimensions with Arrwwid number 

3 nd 

4 ' Z • 

(A Daun tiling lifted into d dimensions has Arrwwid number 3 • 2 d ~ 2 .) 
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Figure 19: A query range for which straightforward proofs lead to overly pessimistic bounds on 
the cover ratio of Dekking's curve. 



7 Discussion and conclusions 

Arrwwid-optimal curves and tilings. This paper shows that in two dimensions the lowest 
possible Arrwwid number is achieved by certain uniform square space-filling curves (such as 
the AR 2 W 2 curve, the Kochel curve, and Dekking's curve) and by a certain rectangular tiling. 
However, for d > 3, no uniform cube space-filling curve can have an Arrwwid number as low 
as the best known rectangular space-filling curve, and the best known rectangular space-filling 
curve does not have an Arrwwid number as low as the best known space-filling curve on a fractal 
tiling. A three-dimensional rectangular recursive tiling with Arrwwid number four or five might 
exist, but so far, it was not found. 

Our lower bounds for space-filling curves only apply to curves based on tilings with simply 
connected tiles (in two dimensions) or with convex tiles (in three dimensions). Lifting the 
restrictions on the shapes of the tiles in these bounds remains a topic for further research 
and I am still pursuing some new ideas on the matter. It may be interesting to investigate the 
relation between the Arrwwid number of a space-filling curve and its maximum multiplicity (the 
maximum number of times any point is visited by the curve). Although the Arrwwid number 
and the multiplicity of a curve may differ (Lebesgue's curve has multiplicity five and Arrwwid 
number four), it is possible that ideas related to the multiplicity of space-filling curves [Uj prove 
useful in deriving bounds on Arrwwid numbers, at least for certain classes of curves. 



Arrwwid numbers versus cover ratios. In principle all proofs of Arrwwid numbers come 
with an upper bound on the cover ratio c. Looking at the proofs in this paper one might 
get the impression that low Arrwwid numbers tend to come at the expense of extremely high 
cover ratios. This may be a shortcoming of the proofs rather than the space-filling curves. 
For example, consider Dekking's curve, which is based on recursively subdividing squares into 
25 squares. Figure [19] shows a disk Q that sticks out just a little bit from the tiles B and E 



with width w, also intersecting A, C, D and F. The approach from the proof in Section 3.3 



would be to cover Q with the parents of these tiles. This results in a cover ratio bound of 
100u) 2 /^tf 2 = 400/7T. However, taking the scanning order into account (see Figure 12), we see 
that Dekking's curve has the property that if B and C are not adjacent in the order, then E and 
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Figure 20: Coil order. 



F must be adjacent, and vice versa. Furthermore, a corner tile and its neighbour are never far 
apart in Dekking's scanning order: there are at most three tiles between them. Note that among 
A, B, D and E there must also be pair that is consecutive in the scanning order. Therefore 
Q can be covered with three fragments containing the six tiles A, B, C, D, E, and F, and at 
most three additional tiles. This results in a cover ratio of at most only 9w 2 /jW 2 = 36/tt (for 
this particular choice of Q). A more detailed case analysis that takes the scanning order into 
account may thus give much better bounds on the cover ratios of large tilings than the bounds 
presented in this paper. 



Worst-case versus average case. Arrwwid numbers only consider the worst-case number 
of tiles or curve fragments that are needed to cover a query range. Therefore lower Arrwwid 
numbers do not necessarily give better disk efficiency on average. It is possible that optimising 
the worst case actually has an adverse effect on the average-case performance. In the case of 
square curves of size four there is some intuition to support this concern: getting an Arrwwid 
number of three requires the use of diagonal connections in this case. For any tile A, let X(A) 
be the set of query ranges Q such that Q intersects A, but not the next tile of the same size in 
the scanning order. The total of |X(A)| over all tiles A gives an indication of the probability 
that after scanning a tile, we must either move the disk head to the next tile that intersects the 
query range, or we scan a tile with false answers only. Since X(A) is larger when A is connected 
diagonally to the next tile than when it is connected orthogonally, diagonal connections may 
lead to decreased performance on average. 



Theory versus practice. In two dimensions, most widely known curves have Arrwwid num- 
ber four (Peano's curve [IT], Hilbert's curve [8], Lebesgue/Z-order [12], Sierpinski/Knopp/H- 
order [16} 118]). while Arrwwid-optimal curves have Arrwwid number three (the AR 2 W 2 -cuvve, 
the Kochel curve, and Dekking's curve, all discussed in Section 3.3). One may wonder if there 
are any practical settings in which it matters much whether three of four fragments are used 
to cover a query range in the worst case. Average-case performance seems to be more relevant 
when comparing these curves: how much time does it cost to scan a query range? I did some 
rough preliminary experiments for certain ratios between seek time (the cost of "jumping" over 
a curve fragment outside the query range) and scanning time (the cost of scanning a fragment, 
relative to its length), but I did not find big differences between the curves. Performance var- 
ied by at most 10%, with the simple coil order (Figure 20) giving the best performance, and 
no apparent advantage for curves with Arrwwid number three. Proper experiments would be 
needed to verify these observations and assess their validity. 

In higher dimensions, the stakes are higher. The recursified shifted cube tiling might seem 
complicated to put to practical use, but one might be able to determine in which tile a point 
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lies by analysing the digits of the coordinates in base-5 notation. This approach may seem 
interesting to explore because in higher dimensions, there is an exponentially increasing gap in 
Arrwwid number between regular hypercube tilings and recursified shifted hypercube tilings. 
On the other hand there is probably also an exponentially increasing gap in cover ratio, in 
favour of the regular hypercube tilings. Further analysis and possibly experiments would be 
needed to make a proper comparison. 

It is therefore unclear yet whether the research described in this paper led to anything that 
may be useful in practice. This research was an investigation into where optimisation of the 
Arrwwid number leads us: What does it take to optimise the Arrwwid number? How limiting 
is it to consider only regular square or cube curves? What are tilings and curves with low 
Arrwwid numbers like? I hope that this work opened some new perspectives on how recursive 
tilings and space-filling curves can be constructed. 
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